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The resummation of multiple soft gluon emission affects the production rate and kinematic distri- 
butions of W + h (where ft is a Higgs boson) and tb pairs at the Tevatron with \/s = 2 TeV. Using 
the Collins-Soper-Sterman resummation formalism, the production rate is enhanced over the next- 
to-leading-order (NLO) prediction by 2-3% for the W + h process, for Higgs boson masses between 
80-120 GeV, and over 3% for the tb process for mt = 175 GeV. After resummation, the tb rate 
changes by 12-13% when mt is varied by ±5 GeV. Various kinematic distributions are presented for 
the individual final state particles and for the pair. The explicit radiation of hard gluons in NLO 
QCD is included also for the tb final state. 
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I. INTRODUCTION 



The Tevatron collider at Fermilab, running at a center- 
of-mass energy y/s — 2.0 TeV with a higher luminos- 
ity, will produce a sizeable number of electroweak gauge 
bosons at large invariant masses, which can subsequently 
decay to top quarks t or Higgs bosons h. Studying the 
properties of t and h are necessary to complete our un- 
derstanding of the Standard Model (SM). The process 
pp — > — > W^h is a promising avenue for detecting 
a light Higgs boson, which subsequently decays h(—* 66) 
HI, or constraining the Higgs boson mass Mh- The pro- 
cess pp — > W + * — > tb (W~* — > tb) has been proposed to 
bound the CKM matrix element Vtb in a kinematic region 
where the parton distribution functions are well known 
0. After the decay t -> bW + (t -> bW~), the second 
process is also a background to the first process, so it is 
interesting to study both in parallel. To determine the 
sensitivity of the Tevatron to Mh and Vtb or to observe 
any new physics effect in the W ± bb system, one has to 
know first the SM prediction for the production rate and 
the kinematic distributions of the final state system. The 
purpose of this study is to quantify the accuracy of the 
SM prediction for the next run at the Tevatron. 

To test the SM or probe new physics, it is important to 
study the kinematics of the final state particles, like h and 
t. An accurate prediction of the rate and kinematics of 
particle production at hadron colliders must include ini- 
tial and possibly final state radiation of gluons in QCD. 
Electroweak corrections are typically much smaller than 
higher order QCD corrections H. The effects of multiple 
soft gluon resummation on the kinematics of on shell W 
bosons and the leptonic decay products have been stud- 



ied in some detail j| . The situation is not so complete for 
off shell boson production at large invariant mass. 
For the processes of interest, the initial and final (if any) 
state QCD corrections at NLO factorize separately. The 
predicted NLO rate for W ± h, with only initial state cor- 
rections, differs from the leading order (LO) rate by 36% 
for Mh =80 GeV. The rate for tb production, with initial 
and final state corrections, differs by 60% for m t =175 
GeV. Initial state soft gluon resummation may affect the 
size of these corrections and modify the kinematic distri- 
butions of the final state as expected from a NLO calcula- 
tion. Furthermore, to accurately predict the kinematics 
of the tb final state at NLO, explicit hard gluon radiation 
must be included. 

Sec.|| contains a review of the Collins-Soper-Sterman 
(CSS) resummation formalism used in this study. We 
closely follow the notation used in Ref. g. In Sec. Ill, 
we present our numerical results for qq' — ► W + * —> W^h 
production using this formalism. We show how to incor- 



porate the effects of final state radiation in Sec. |V, and 
present numerical results for qq 1 — » W + — > tb. Finally, 
Sec. [v| contains our conclusions. 



II. THE CSS RESUMMATION FORMALISM 



Soft gluon resummation has been applied successfully 
to predict the rate and kinematics of electroweak gauge 
boson production at hadron colliders For W ± 

bosons, the CSS formalism has been applied only to 
leptonic final states. The same formalism can be di- 
rectly applied to W + — * W + h production. However, 
for W + — > tb production, there are also singularities as- 
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sociated with the final state radiation which need to be 
handled separately. 

The resummed differential cross section in the CSS for- 
malism for the production of gauge boson V at a hadron 
collider with a center-of-mass energy is 



/ dajfufo -> V(-> did 2 )A") \ 
\ dQ 2 dy dQ\, d<pv d cos 9 d<fr J re 







d 2 b 



JQT-b 



967r 2 s(Q 2 -M 2 ) 2 + M 2 T 2 {J (2tt) 2 

x w jk{b*,Q,x 1 ,x 2 ,9,(j), C 1 ,C 2 ,C 3 ) 

xF™ > (b,Q,x 1 ,x 2 ) + Y{Q T ,Q,x 1 ,x 2 ,6,cf>,C i )^, (2.1) 

where Q, M, r, y, Qt and <py are the mass, on shell mass, 
width, rapidity, transverse momentum, and azimuthal 
angle of V in the laboratory frame, 8 and <f> are the 
polar and azimuthal angle of the final state particle 
d\ [see Eq. fl2.8| )] in the Collins-Soper frame |9j, x\ = 
(Q/y^) exp(y]~z 2 = (Q/y/s)exp(-y), and (3 = 2p*/Q, 
where p* is the momentum of di and d 2 in the V rest 
frame. The renormalization group invariant Wjk is given 
by 

W jk (b,Q,x 1 ,X2,8,<f>,C 1 ,C 2 ,C 3 ) = 
exp{-S(&,Q,Ci,C 2 )}|^ fe | 2 

x <^ [(Cji <g> fi /hl ) {x x ) (Cfem ® fm/hi) (^2) 

+ {Ckl ® fl/hi) {Cjm ® /m/Tia) O^)] 

X 9lfp iS 3 k(Q^ COs0 ^ rn di,m d2 ) 

+ [(Qi ® /j/fci) fal) (Cfcm ® f m /h 2 ) (X 2 ) 
- {Ckl <8> /i/ftj (iCl) (Cjm <g> / m /fe 2 ) (x 2 )] 

X 5i/iW/ fe (Q,cos6»,m dl ,m d2 )|, (2.2) 

where Vjk is a CKM matrix element for the initial state 
partons j and k, S is a Sudakov form factor, m^m^ 
are the final state masses, T-L s ^ is a symmetric (anti- 
symmetric) function in cos 0, g\ = fs(A) 
is a coupling for the symmetric (antisymmetric) func- 
tion, and <g> denotes the convolution integral of the Wil- 
son coefficients C and the parton distribution functions 
(PDF's) 
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{Cji®fi/ hl ) [xi) = 

^CAfA^^CM) 

i si si o 

*fi/hi(£utJ>= -r)- 



(2.3) 



For the processes under consideration, j — q,k — q , 
excluding t quarks. The dummy indices I and m sum 



over quarks and anti-quarks or gluons, and summation 
on double indices is implied. The angu lar f unction TL S ^ 
and coupling constants fs{A) m Eq. fl2.2j ) for the di — 
W + , d 2 = h final state are 

(M^/Q 2 ) [2 + (p*/M w ) 2 (l - cos 2 9)] , (S) 

0, J (A) 

f 2 = V2G F M 2 v J 2 A = 0, 

and for the d\ = t, d 2 = b final state are 

3 [(1 - (m 2 - m 2 ) 2 /Q 4 + (2p*/Q) 2 cos 2 8)] , (5) 
3[(4p7Q)cos0], (A) 

fl = f 2 A = al\v tb \ 2 . 



The Sudakov form factor S(b,Q,C u C 2 ) in Eq. fl2_l|) is 
defined as 



c- 



d/2 2 



In 



/i 2 



A(a s (p,)) +B(a a (p)) 



(2.4) 



where a s is the strong coupling constant, and the func- 
tions A, B and the Wilson coefficients Cji were given in 
Ref. jsj. The constants Ci, C 2 and C3 = /16 arise when 
solving the renormalization group equation for the Wjk- 
The canonical choice of these renormalization constants 
is C\ = C3 = 2exp(— 7b) = 60 and C 2 = 1, where 
is the Euler constant. In general, the choice C\ = C 2 bo 
and C3 = bo eliminates large constant factors in the ex- 
pressions for the A, B and Cjk functions. As shown in 
Eq. (2.4), the upper limit of the integral for calculating 
the Sudakov factor is p, = C 2 Q, which sets the scale of 
the hard scattering process when evaluating the renor- 
malization group invariant quantity Wjk, as defined in 
Eq. (2.2). The lower limit \x = C\/b = bo/b determines 
the onset of nonperturbative physics. When integrated 
over the impact parameter b, the W term is referred to 
as the CSS piece. 

The Y piece in Eq. 12. 1|) is defined as 



Y(Q T ,Q,xi,x 2 ,6 1 4>,C i ) 



1 d£i 
x si 



1 J<- 00 
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IT 
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(2.5) 



where the functions iJ.„ only contain contributions less 
singular than Q^, 2 x (1 or ln(Q 2 / Q 2 *-)) as — > in the 
fixed order perturbative calculation. The perturbative 
expansion is optimized by choosing C4 = 1, so that the 
scale of the Y piece is Q. 

In Eq. (2.1), the impact parameter b is to be integrated 
from to 00. However, for b > 6 ma x) which corresponds 
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to an energy scale much less than l/b max , the QCD cou- 
pling a s becomes so large that a perturbative calculation 
is no longer reliable. The nonperturbative physics in this 
region is described by an empirically fit function F NP 
with the general structure 



exp < — In 



F J N ) r(b,Q,Qo,x 1 ,x 2 ) 
-h k /h 2 (x 2 ,b) 



(2.6) 



The functions hi, hj/h x and hp-i^ cannot be calculated 
using perturbation theory and must be measured exper- 
imentally. Furthermore, W is evaluated at 6*, instead of 
6, with 



6* = 



y/l + (ty&ma*) 2 



(2.7) 



such that 6* never exceeds & ma x- 

To obtain the kinematics of the final state products d\ 
and d 2 from V* — ► did®, we transform the four momen- 
tum of di (= p^) and d 2 (= p M ) from the Collins-Soper 
frame to the laboratory frame. The resulting expressions 
are 

" " % ("77 + sin 6 cos $ Xil + sin0 sin «/>y + cos eZ' 1 ), 



P' 



2 V Q 



q M = [Mt cosh y, Qt cos 4>v , Qt sin (f>y , Mt sinh y) , 
X" = 



Mt 



(2.8) 



with q» = {q°, q\q 2 ,q% Q± = {q° ± q 3 )/V2, Q = 
y/f, M T = VO 2 + Or, V = ln(Q+/Q_)/2, rf = 
(1,0,0, 1)/V2, n v = (l,0,0,-l)/\/2, and e 0123 = -1. 
In the Collins-Soper frame of reference, the W* is at 
rest and the z axis is defined as the bisector of the an- 
gle formed by p% x and — j>h 2 , where p^ and ph 2 are the 
momentum of the initial state hadrons. As suggested by 
the notation, X M , Y^, are the transformations of the 
x— , y—, z— components of d%, respectively, from the CSS 
to the laboratory frame. Similar expressions can be de- 
rived when there are more particles in the final state as 
discussed in Sec. IV 



In this section, we present numerical results for the 
qq' — > W + * — > W + h process after applying the resum- 
mation formalism outlined in the previous section at the 
upgraded Tevatron with \fS = 2.0 TeV. For these re- 
sults, we have assumed nit — 175 GeV and m& = 5GeV 
and several values of the Higgs boson mass Mh- In the 
following, distributions will be shown for W + h produc- 
tion only. This allows us to exhibit asymmetries along 
the beam axis. The distributions for W~h production 
can be obtained by reflection about rapidity y = 0. 

As explained in the previous section, the CSS piece de- 
pends on the renormalization constants C\,C 2 — C\/bo 
and C3 = bo, as well as a few other implicit parameters. 
The choice of C 2 indicates that the hard scale of the pro- 
cess is C 2 Q- For this study, we use C 2 — 1, which fixes 
the scale at the off shell gauge boson mass. We also use 
the CTEQ3M NLO PDF's §J, the NLO expression for 
a s , and the non-perturbative function Ji"l| 

F^{b,Q,Q , Xl ,x 2 ) = 
1 2 i2 i_ / Q 



exp <^ - gib - g 2 b In 



2Qo 



-gig?,b In (IQQxix 2i 



(3.1) 



where gi = 0.11 GeV 2 , g 2 = 0.58 GeV 2 , g 3 = 
-1.5 GeV" 1 and Q = 1.6 GeV. The choice 6 max = 
0.5 GeV" 1 was used in this fit. Finally, the CSS piece 
is fixed by specifying the order in a s of the A, B and 
Cjk functions. We adopt the notation (M, N) to repre- 
sent the order in a s of ^,5^ and cjf } . The choice 
(1, 0), for example, means that A and B are calculated to 
order a s , while Cjk(z) is either or 5(1 — z) depending 
on j and k. 

The CSS piece alone gives an accurate description of 
the kinematics of the off shell W + boson for Qt <C Q- 
However, because of the soft gluon approximations made 
in the resummation, the CSS piece becomes negative 
when Qt ~ Q- The regular Y piece accounts for terms 
neglected in W at any fixed order in perturbation theory. 
The sum CSS + Y can become negative for Qt ~ Q, but, 
when integrated over Qt from to Q, returns the NLO 
result for the total cross section and Q distribution. We 
have checked this agreement numerically. Since the Qt 
distribution is a physical observable, negative differential 
cross sections are not allowed. Therefore, when CSS + 
Y no longer gives an accurate description of the physics, 
one should use a prediction at a fixed order in a s , since 
the perturbation series in a s gives a better convergence 
for Qt ~ Q- Switching between the two formalisms re- 
quires a matching prescription. Our prescription is to 



III. THE NUMERICAL RESULTS OF INITIAL 
STATE RESUMMATION 



* These values were fit for CTEQ2M PDF and C 2 =l, and 
in principle should be refit for CTEQ3M PDF. 
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switch smoothly to the fixed order differential cross sec- 
tion once CSS+y becomes smaller than the fixed order 
prediction as Qt — Q- The cross section after matching, 
then, gives a different prediction than the NLO result. 
We refer to a full resummed calculation, including the Y 
piece and matching, as CSS(M,N), where (M,N) specify 
the order of the A, B, and C functions. 

Part of the next-to-next -to-leader-order (NNLO) 
prediction for W* prod uctio n can be included through 
and in Eq. (Q §. When matching to 

the CSS(1,1) (CSS(2,1)) calculation, we use the NLO 
(NNLO) perturbative prediction for high transverse mo- 
mentum boson production |l2| . We apply the results 
of our calculations to study the total rate and the kine- 
matic distributions of the W + h pair and the individual 
W + or h produced in hadron collisions. These do not 
include any final state effects. 

The effects of resummation on the total production 
rate are illustrated in Table |[ which compares LO, NLO, 
the CSS formalism with A, B, and C functions of or- 
der (1,1) and matching, and the CSS formalism of order 
(2,1) and matching. The LO result uses the CTEQ3L 
PDF's. The CSS(1,1) and CSS(2,1) predictions are in- 
distinguishable once matching is used, and reveal a slight 
enhancement over the NLO prediction. We interpret the 
difference between the NLO and CSS predictions as a 
gauge of the intrinsic theoretically uncertainty in the to- 
tal rate up to order a 2 . Referring to Table |, higher order 
corrections to the NLO prediction could be as large as 2- 
3%. There will also be some uncertainty from varying 
the PDF, but we have not investigated this. 

The effect of resummation on various kinematic distri- 
butions are displayed in Figs. [l]-^[ Distributions for the 
individual W + and h are basically changed in rate, but 
not in shape, compared to the LO and NLO predictions 
(see Figs. |]-|]). On the scale of the figures, there is not 
a significant difference between CSS(2,1) and CSS(1,1), 
so we do not include CSS (1,1) in these figures. However, 
those observables which depend on the kinematics of the 
pair are sensitive to the order of resummation. For the 
purposes of generating NLO distributions, a Qt cut of 
.8 GeV is imposed to separate the real emission from the 
virtual corrections. Qt is the transverse momentum of 
the virtual W + boson, or, equivalently, the W + h pair, 
since we do not include any final state corrections. Near 
this cut, we have little confidence in the ability of the 
NLO calculation to accurately describe the Qt distri- 
bution. Even far from this cut, however, the NLO Qt 
distribution grossly underestimates the Qt distribution 
as expected from resummation (see Fig. [y). At high Qt, 
CSS(2,1) predicts an enhancement over CSS(1,1), which 
means more hard gluon radiation. Also, CSS(2,1) pre- 
dicts a large azimuthal angle separation (A<fi w h ) be- 
tween W + and h (see Fig. ||). 



IV. INITIAL STATE RESUMMATION WITH 
EXPLICIT FINAL STATE RADIATION 



Using the same approach discussed in Section III 



we 

have also studied the effects of soft gluon resummation 
on the process qq' — > W + * — > tb for several values of 
m t — 170, 175, 180 GeV. The presence of hadronic fi- 
nal states, however, adds a complication. In Ref. [13[ , 
the CSS formalism was applied to the production of tt 
pairs at the Tevatron. The effects of final state radiation 
were incorporated in this result by resumming all of the 
final state logarithms that behave like Q^ 2 ^{Qt/Q) 2 
and adding these coherently to those similar terms from 
the initial state. This was justified on the grounds that 
color links the initial state to the final state at NLO. 
Also, previous studies of top quark decay at NLO show 
that most gluon radiation occurs off the b quark line, so 
that final state radiation off the t quark line is a smaller 
effect [Q. For the production of a hadronic final state 
through an s channel W* boson, there is no interference 
at NLO between Feynman diagrams with initial and fi- 
nal state QCD corrections, so that the effect of initial 
state corrections can be calculated separately from final 
state corrections. To calculate the total cross section for 
W* — > tb, we resum the initial state multiple soft gluon 
radiation and include the effect of QCD radiation at NLO 
for hadronic final states by multiplying the production 
rate for pp — > W + * — > tb by a factor of Kf sr (Q 2 ), a func- 
tion of only Q 2 once the masses of t and b are fixed [jl5| . 

The effects of initial state resummation with the fi- 
nal state QCD correction K{ sr on the total production 
rate are summarized in Table ||, which compares the 
LO, NLO, CSS(1,1), and CSS(2,1) predictions for var- 
ious values of m t . For m t = 175 GeV, the enhancement 
from the CSS formalism over NLO is about 3.5%, which 
is an estimate of the effects of higher order corrections. 
The enhancement over the LO result is 66%. As a com- 
parison, the variation in rate from a change in nit by ± 
5 GeV around a nominal value of 175 GeV is 12-13%. 
Our NLO correction from initial state radiation, which is 
37% for m t — 175 GeV, agrees with Ref. ||. However, 
our NLO correction from final state radiation, which is an 
additional 17% enhancement, is about 3.4% higher than 
Ref. H. The enhancement of the resummed result over 
the NLO result is based on our calculation of the NLO 
rate. A study of the kinematic distributions of the final 
state particles t and b at NLO in QCD requires a sepa- 
ration of the states with and without explicit hard gluon 
radiation. While Kf ST can predict the change in rate from 
final state radiation, it has integrated out the kinematics 
of the final state particles. To remedy this situation, so 
that we can study the explicit final states tb and tbg, we 



calculate the correction function K^{Q 2 ,E\ 



f m ) for the 



exclusive process W + — > tbg using a minimum gluon 
energy E™ ln as a regulator. The effect of the regulator 

on the exclusive process W + * — > tb can be calculated us- 



mgK$(Q 



2 ^min 



) = k {si (Q 2 )-k£ ) (q 2 ,e ! 



™ m ). This 
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is a practical method as long as the value of P™ ln trans- 
formed to the laboratory frame does not exceed the low- 
est jet energy used by an experiment in an analysis. In 
our numerical results, P™ m = 1 GeV. fj] 

Figs. p|-|rT|show the effects of initial state resummation 
and final state NLO QCD radiation on the kinematic dis- 
tributions of the t, b, and the pair tb. The Qt distribution 
in Fig. |^ is still the transverse momentum of the virtual 
W + * boson, so it does not depend on the separation of 
final states with and without hard gluon radiation. For 
this study, the t and b are produced on shell and treated 
at the parton level. No jet definition is used to combine 
the hard gluon with t or b partons. We present results 
for the resummation calculations CSS(2,1). On the scale 
of these figures, there is no significant difference between 
CSS(2,1) and CSS(1,1). For comparison, we show the 
distribution expected using two body kinematics for tb 
but still using the full final state correction Kf SI . The 
difference between this distribution, denoted as 2-body 
in the figures, and the resummed curve indicates the size 
of the "hard" gluon effects, that is, the contribution from 
tbg for a fixed P™ m . When we feel a comparison is useful, 
we also include the LO prediction. For the kinematics of 
the individual t or b, the effect on the shape of distribu- 
tions from final state radiation is important, more so than 
from the resummation (see Figs. |7|-|lT|), and we do not see 
a noticeable difference between the NLO and resummed 
curves with our choice of P™ m . The kinematics of the 
tb pair, however, are sensitive to the initial state gluon 
radiation (see Figs. |u] and |l2]) . This can be understood as 
follows. Similar to the W^~h final state, intial state gluon 
radiation only affected distributions of the tb pair. Why 
is this so? This is mainly determined by the kinematics 
and the nature of the radiation divergence for the ini- 
tial and final states. In the initial state, the partons are 
collinear with the beam axis, and the leading divergence 
is both soft and collinear. As a result, the Qt distri- 
bution (see Fig. ||) peaks around a few GeV, while the 
rapidity of the pair is very central. The LO distribution 
for P T , in contrast, peaks near 30 GeV. Therefore, there 
is little effect of the Qt smearing on P T . The same ar- 
gument can be applied to the t rapidity distribution y l . 
The final state radiation, on the other hand, only pro- 
duces soft divergences, and the final state particles tend 
to be central and not collinear with the beam. Therefore, 
distributions with respect to the beam axis, such as P T 
and P T , are more sensitive to the final state effects. 



' As discussed in Ref. |J, if the kinematic cuts imposed 
to select the signal events are correlated with the transverse 
momentum Qt of the virtual W boson, then the distributions 
predicted by the NLO and the resummed calculations will be 
different. The difference is the largest near the boundary of 
the phase space, such as in the large rapidity region. 



V. CONCLUSIONS 



We have studied the effects of QCD resummation in the 
CSS formalism on the production of W + * bosons at the 
Tevatron and the kinematics of the final state products 
W + h and tb. Understanding the W + h production rate 
is important for quantifying the Higgs boson discovery 
potential of the Tevatron. We find an increase in rate 
over the LO result of 40% for M h = 80 GeV, and 29% 
for Mh = 120 GeV. This is only slightly higher than the 
NLO prediction, even though we have included part of 
the NNLO correction. Therefore, we estimate the effect 
of higher order QCD corrections to be 2-3% over the 
NLO prediction. The shape of kinematic distributions 
of the W + and h are not significantly altered from the 
LO or NLO predictions. The kinematics of the W + h 
pair, however, can only be accurately predicted using a 
resummed calculation, and might be useful to distinguish 
signal from background. 

Understanding the tb production rate is important for 
the determination of the CKM matrix element Vth- We 
have included the explicit radiation of final state gluons 
at NLO, so that we predict the kinematics of tb and tbg 
events. We find an increase in rate over the LO result 
of 66% for m t = 175 GeV, and we estimate the effect 
of higher order corrections to be more than 3% over the 
NLO result. The total production rate changes 12-13% 
when nit is varied by ±5 GeV. The effects of final state 
radiation are more important than initial state resum- 
mation for the kinematics of the individual t or b. With 
explicit gluon radiation, the Pr spectrum of both t and 
b are harder, so that the b tagging efficiency should be 
higher. The kinematics of the tb pair are also interesting 
observables. 
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FIG. 1. The Q T spectrum for W+h Events 







Higgs Mas 


s (GeV) 




Order 


80 


100 


120 


CSS (2,1) 


.305 


.154 




.083 


CSS (1,1) 


.304 


.154 




.083 


NLO 


.296 


.149 




.081 


LO 


.218 


.109 




.059 



TABLE I. Total cross section in pb for W + h production at 
the Tevatron in various approximations for M h =80,100, and 
120 GeV. 



Order 


Top Mass (GeV) 
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180 


CSS (2,1) 


.541 


.478 
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CSS (1,1) 


.541 


.477 




.423 


NLO 
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.462 




.409 


LO 


.327 


.288 




.255 



TABLE II. Total cross section in pb for tb production at 
the Tevatron in various approximations for mt =170, 175, and 
180 GeV. 
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FIG. 2. The P% spectrum for W + h Events 
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FIG. 4. The Ay w h spectrum for W + h Events FIG. 6. The Q T spectrum for tb(g) Events 
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FIG. 11. The cos6>*^ spectrum for tb(g) Events 
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FIG. 12. The A(f> tb spectrum for tb(g) Events 



0.012 



M t =175 GeV 




20 40 60 80 100 120 140 160 180 200 

Pj (GeV) 



M t =175 GeV 
CSS(2,1) 
CSS(1,1) 

LO 
2 Body 



